Copyright © 1981 American Telephone and Telegraph Company 

The Bell System Technical Journal 

Vol. 60, No. 10, December 1981 

Printed in U.S.A. 



Study of a Time-Compression Technique for 

TV Transmission Using a Chirp Filter and 

Envelope Detection 

By K. Y. ENG and B. G. HASKELL 

(Manuscript received January 19, 1981) 

We study a time- compression (or expansion) technique for possible 
application in communication signal processing, e.g., broadcast- 
quality tv transmission through satellites. The method uses a linear 
chirp, a linear dispersive filter realized by surface acoustic wave 
devices and an envelope detector. This technique is heuristic and can 
be viewed as a quasistationary model of the fm wave involved. 
Numerical results show that excessive distortion is created, and its 
application to tv transmission is not suitable unless some kind of 
equalization is provided. One such form of equalization is the chirp 
transform processor which involves considerably more complexity. 
Simpler equalizations may be possible but do not seem to be straight- 
forward. 

I. INTRODUCTION 

We study a time-compression technique motivated by the long- 
standing interest in transmitting multiple broadcast-quality color TV 
signals through a single satellite transponder, i.e., a usable rf band- 
width of 36 MHz in a communications satellite such as COMSTAR. 
This can be done by the use of frequency division multiplexing (fdm). 
However, the nonlinearity of the transponder can cause serious intel- 
ligible crosstalk and intermodulation interference between the fm 
carriers unless the satellite power amplifier is backed off substantially. 
Such a backoff, in turn, leads to a reduction in the downlink carrier-to- 
noise ratio. As a result, there exists an optimum trade-off between the 
crosstalk and s/n's, which limits the overall system performance, and 
achieving broadcast-quality TV transmission becomes difficult. 

It is possible to time compress each scan line of a color tv signal by 
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the use of a linear chirp, a linear dispersive filter (ldf) and an envelope 
detector. Two or more time-compressed scan lines from different, but 
synchronized, TV signals can then be time multiplexed together in the 
time duration of an ordinary TV scan line. This concept of time- 
compression multiplexing (tcm) is not new, 1,2 but recent advances in 
fast analog-to-digital converters, digital-to-analog converters and 
charge-coupled devices have greatly facilitated the implementation of 
time compression or expansion. However, because of their present 
limitations on bandwidth and speed, time-compression factors for 
achievable TV signals are only 2 or 3. For large time compressions, 
ldfs realizable by surface acoustic wave (saw) technology are prom- 
ising candidates because of their large time-bandwidth property. In 
addition to high speed, the TV application also requires extremely high 
signal fidelity. There are other applications on the high-speed time 
expansion (or compression) of waveforms where the distortion require- 
ment is less stringent than the tv transmission case. In the specific 
case of multiple tv transmissions through a single satellite transponder, 
there are many advantages in using tcm, e.g., higher transponder 
efficiency, no intermodulation, no crosstalk, possible compatibility 
with time division multiplex (tdm) operations, etc. The crucial ques- 
tion is, of course, how much distortion the compression/expansion 
process would introduce on the signals. This paper gives both analysis 
and numerical examples that illustrate the method. 

The study revealed that considerable distortion is introduced by 
these operations, and its application to broadcast-quality tv transmis- 
sion would require saw filter performance beyond the present state of 
the art. However, if the distortion requirement can be relaxed some- 
what, then the present approach is advantageous because of its sim- 
plicity. On the other hand, if high signal fidelity is required, then some 
kind of equalization is needed for the present technique. This has 
motivated the study of an extension of the present method which is 
capable of producing high signal quality with saw filter requirements 
within the present state of the art even at compression factors of 10 or 
more, but at the expense of higher complexity. This latter development 
is not discussed here but is covered in Ref. 3. The remainder of the 
paper covers the theoretical analysis and computer simulation. How- 
ever, the discussion of either subject by itself is not adequate for the 
complete understanding of the system. The theoretical analysis estab- 
lishes that although the basic concept was derived heuristically 
through physical interpretations, it can be viewed as a quasi-stationary 
approximation to the time-compression process. The computer simu- 
lation, on the other hand, provides the quantitative results that lead 
to the conclusion that the resulting distortion is excessive for today's 
saw filter parameters. 
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II. THEORETICAL ANALYSIS 

In this section, we describe and analyze the proposed compression 
method using tv as an example. We first describe a heuristic argument 
of how the technique is supposed to work. We then derive the impulse 
response of a general ldf — the understanding of which is important to 
the subsequent analysis and simulation. A brief step-by-step analysis 
of the compression process is shown, and its result reveals that the 
technique can be interpreted as a quasi-stationary approximation of 
the chirp signal. The mathematical expressions describing the time 
compressor are complicated; thus, numerical results are obtained using 
a computer simulation discussed in Section III. Various other proper- 
ties are also discussed. 

2.1 A physical interpretation 

A block diagram is shown in Fig. 1. The input signal v(t) consists of 
successive scan lines, each with a duration of Ti seconds, and the 
voltage is biased to be positive. It is multiplied synchronously by a 
periodic chirp signal c(t) with a center frequency f and a chirp range 
of A/ c . The instantaneous frequency of c(t) sweeps linearly from 
(/o - A/c/2) to (fo + A/c/2) over each scan line duration Ti. The lowest 
frequency of the chirp signal is assumed to be much greater than the 
highest frequency in the tv signal. The input x(t) to the ldf is then 
an amplitude-modulated chirp waveform. For simplicity, let us restrict 
our attention to a single scan line, say (0, TV). In this interval, x{t) 
chirps from {f - A/ c /2) to (f + A/" c /2) with the tv signal as the 
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Fig. 1 — Time-compression filter. 
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envelope modulation. As for the ldf, we assume that it has a band- 
width A/ centered at f , and A/> A/" c . 

Again, for the purpose of a simple illustration, let us assume that 
A/ = A/ c and, over its passband, the ldf has a constant gain and a 
linear group delay characteristic decreasing from At to zero, where At 
is the delay dispersion of the ldf. At the time instant t = 0, x(t) has 
an instantaneous frequency (/ - A/ c /2) and an envelope magnitude 
proportional to v(0). This "envelope piece" is delayed by At as it 
transits through the ldf. Similarly, at t = Ti, the instantaneous 
frequency of x(t) is the highest chirp frequency which gives a zero 
delay for the passage through the ldf. In between these two end 
points, the envelope of x(t) is delayed linearly from At to zero. 
Equivalently, the envelope of x(t) over the interval (0, T t ) is time 
compressed into (At, Ti) in the ldf output. An envelope detector is 
then used to retrieve the time-compressed TV signal. Similar compres- 
sion occurs for all the other scan lines, and as a result, the envelope 
detector output consists of a sequence of compressed TV bursts. If we 
denote the duration of each burst by T c , the ratio (Ti/T c ) is called the 
time-compression ratio (tcr). 

It is easy to show that, in general, for a given set of Ti, At, A/", and 

A/c, 

Tc = T/ ~a7 a/c ' (1) 

where < T c < Ti, A/ c < A/", and the compression ratio is 

It is clear from the above description that time expansion is also 
possible by the use of an increasing delay characteristic in the ldf, 
and in such a case, the filter will become a time-expansion filter. 

2.2 Impulse response of a general LDF 

We consider the general case of an idealized ldf as shown in Fig. 2a. 
The bandwidth of the filter is A/ 7 centered at f . The group delay varies 
linearly over the passband as shown in the diagram. The transfer 
function of the filter, using analytic-signal notation, is 



H(f) = 



exJ-flw (f-fo)r + if J 0) +fr J, 



-¥-'-*«¥■ 



0, elsewhere, (3) 
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Fig. 2— Characteristics and impulse response of an ldf. (a) Transfer function, 
(b) Impulse response. 



where 



A/ 

At' 



(4) 



to is the group delay at f , and <f>i is a constant. The inverse Fourier 
transform of H(f) gives the impulse response h(t). The derivation for 
h(t) is given in Appendix A. Neglecting some small envelope and phase 
perturbations, a good approximation (sufficient for our present consid- 
eration) for h(t) is 



h(t)=- 



COS 27T 



r AA a , 
/0 "T ' + 2* +** 



At At 

To ~T-'- To + T' 



0, elsewhere, 



(5) 



where the multiplying constant has purposely been dropped, and <fc is 
a constant. A sketch of h(t) is shown in Fig. 2b. It is a linear chirp 
waveform starting at t - At/2 and ending at t + At/2, with the 
instantaneous frequency varying from f - A//2 to f + A//2 at a chirp 
rate of a. 

2.3 Analysis of time compression 

In this analysis, we again restrict our discussion to a single scan line 
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for simplicity. We denote this input scan line by A(t) in the interval 
(0, Ti). Referring to Fig. 1, the linear chirp signal is given by 



c(t) = cos 2tt 



fo- 



H 



P < 



+ <t>3 



0<t<Ti, 



where p is the chirp rate defined by 



'-$ 



(6) 



(7) 



and 08 is a constant. The parameters f and A/ c are the chirp center 
frequency and deviation, respectively. There are also two implicit 
assumptions: (i) f » A/ c ; and (ii) f » the highest frequency in the tv 

signal. 
The ldf input is the product of the tv input A(t) and c(t), i.e., 



x(t) = A(t)cos 2ir 






0<t<Ti. (8) 



To obtain the ldf output, we can either use the time-domain approach 
by convolving x(t) with the impulse response h(t) of the ldf, or use 
frequency domain analysis by multiplying the Fourier transform of 
x(t) by the ldf transfer function and then performing an inverse 
transform. The former method is much simpler and is presented here. 
The latter is tedious, offers no additional insight, and is, therefore, 
deleted for brevity. 

The ldf is assumed to have an extended passband over (f Q - A//2, 
fo + A//2), where A/"> A/ c . Its delay characteristic is shown in Fig. 3. 
Note that the delay slope is given by 




FREQUENCY (/I 

Fig. 3 — Delay characteristic of the ldf. 
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A/" 



(9) 



Also note in Fig. 3 that the delay at f + At/2 is zero, as some constant 
delay through the device has been dropped for simplicity. Using the 
result of Section 2.2, the impulse response of the ldf is 



h(t)=- 



COS 2?T 



fit --t 2 + <f>4 



, < t < At, 



0, elsewhere, (10) 

where /i A f - A//2, and <> 4 is a constant. The ldf output is then 

y(t) = I x(r)h(t - r)dr 



A(r) 



X^COS277 



®(t) + (at-bf c -f 4 + f 3 )T+(p-a)- 



COS 27T 



$(0 + (2/ -at + f 4 - h)r + (/? + a) - 



rfr, 



< * < T, + At, 

where the limits of integration are defined by 

Tx A max(0, t - At), 
T 2 A min(«, T,), 

/" 3 and U aie defined in Fig. 3, and 

0(0 A fa - ^t 2 + <$>< . 



(11) 

(12) 
(13) 

(14) 



A constant phase term has been dropped in (eq. 11), and we will 
neglect all unimportant constant multipliers and phase shifts in sub- 
sequent discussions. The integral of the second cosine term in eq. (11) 
can be discarded because of the high frequency component 2/" t in the 
integrand. Therefore, y{t) becomes 



y(t) 



■r 



j4(t)C0S 277 



®(t) + (at-bf c -f 4 +f 3 )T 

r 2 ~\ 
+ (P-a)j dT, < t < Ti + At. (15) 
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Let us examine the above integral expression carefully. Using analytic- 
signal notation, we rewrite it as 

y(t) = I A(r) expj j2ir\ $(t) + (at - A/ c - f 4 + f 3 )r 



+ &-«)- 



= expy'2 



ftirQ(t) 



A(t) expj >2t7 



(a* - A£ - U + /s)t 



+ OS - «) g K f 

< t < T, + At. 
In this form, we can define a complex envelop for y(t ) as 

A y (t) A I A(r) expj >2tt 



(16) 






(at-bfc-f4 + h)r 



+ (/? ~ «) J 



rfr k < * < T/ + At. (17) 



In the above, note that A(t) is a slow-moving function while the 
exponential term contains a highly oscillatory chirp given by the 
derivative of the bracketed argument with respect to t, i.e., 



fi(r) = (at - A/c - U + U) + (i8 - «K Ti < t < T 2 



(18) 



This can also be obtained by simply taking the difference of the 
instantaneous frequencies of x(t) and h(t - t) in eq. (11). The convo- 
lution integral involved is illustrated in Fig. 4, where we show x(t) and 
also the corresponding /|(t) at a fixed t = ft. It can be seen that A y (t\) 
is given by an integral over (T x , T 2 ) of a linear chirp waveform at a 
chirp rate of (/? - a) and with an envelope modulation A(t). Further- 
more, the chirp frequency ft inside (T h T 2 ) may vanish at some t as 
shown in Fig. 4. In such a case, the value of y(t\) is dominated by the 
integral eq. (15) over the small interval surrounding that t, where /i 
goes to zero. This is, of course, the well-known quasi-stationary ap- 
proximation. The approximation is good if the chirp rate (/? - at) 
and the interval (7\, T 2 ) are large, and A(r) variation is slow by com- 
parison. Using this approximation, at t — ft inside the valid interval 
(0, T, + At), 



A,(ft) * kA(ri), 



(19) 



2380 THE BELL SYSTEM TECHNICAL JOURNAL, DECEMBER 1 981 




INSTANTANEOUS 
FREQUENCY 



Fig. 4 — Convolution integral. 



where k is a constant, and ti is obtained by solving eq. (18) with fi set 
to zero and t = ti, i.e., 

ah - A/c - f 4 + f 3 



T] = 



P-a 



(20) 



This quasi-stationary approximation is indeed equivalent to the phys- 
ical interpretation described in Section 2.1. As a check, let us derive 
the tcr from the approximation above. We know that A(t\) is nonzero 
only if < ti < T/. The corresponding t } can be solved for using eq. 
(20), and the end points of t\ constitute the interval T c , i.e., 



T c = 



(A/c + U ~ h) ~ TAP - a) (A/ c + f A ~ H) 



a, 



(21) 







Therefore, 

TCR=/l-^l , a>/?, (22) 

which agrees with eq. (2) with the substitutions of a and ft according 
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to the definition eqs. (4) and (7), respectively. See Section 2.4 for a 
continued discussion of A y (t). 

To finish our analysis of the time compression, we return to y(t) in 
eq. (15) and expand the cosine term: 

y (t) = ;y c (*)cos 2nQ(t) - y s (t)sin 2irQ(t), (23) 



where 



•T 2 

y c (t) A I A(t)cos 2tt 



and 



'T 2 

y s (t) A I A(r)sm 2-n 



T 2 ' 

(at - ^fc - U + fz)t + (P - a) - 



T 2 

(at - A/ c - U + fs)r + (0-a)- 



dr (24) 



dr. (25) 



In eq. 25, y c (t) and y 8 {t) can be viewed as the in-phase and quadrature 
components of y(t) after synchronous demodulation. The envelope 
detector output is then 

z(t) = [y 2 c (t)+y 2 s(t)] i/2 . (26) 

2.4 Some fundamental properties 

We have just derived mathematical expressions for the time-com- 
pression process. These expressions are too complicated for easy 
interpretation. However, we have demonstrated that the technique 
works if a quasi-stationary approximation is made for the chirp 
waveform. In other words, the instantaneous frequency of the chirp 
wave could be used as if it were a stationary carrier frequency in a 
steady state analysis. 

Making such an assumption, it is seen from eq. (2) that infinite 
compression, TCR - oo, results if A/ = A f c and At - Ti (a = 0). But 
from eqs. (15), (24), and (25), it is obvious that the ldf output after 
synchronous detection (for a = fi) will actually become the Fourier 
transform of the input envelope A(t). This can also be recognized as 
the well-known chirp transform or real-time transform 4,5 commonly 
used in chirp radar and saw processors. Therefore, the quasi-stationary 
model is invalid for this case. 

A case where the quasi-stationary approximation clearly holds is 
where A/ c -> °° and a~ l -* 0, i.e., the chirp range is very large, and the 
delay slope of the ldf is close to zero. The result is, of course, a very 
slight compression, i.e., the TCR is slightly larger than 1. Therefore, 
without doing any specific calculation, we see that the quasi-stationary 
assumption is valid, at best, for small tcrs, and it breaks down 
somewhere between tcr = 1 and oo. Since our practical applications 
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require tcr > 2, the case of tcr = 2 becomes most interesting and is 
investigated in the simulation later. 

Let us now assume the quasi-stationary model and see what kind of 
distortion would result. The input to the ldf can be expressed as a 
multitone am signal: 

x(t) = I 1 + £ mjcosujt Jcosojc^, (27) 

where uj are the angular modulating frequencies, and w c is the angular 
carrier frequency at some fixed instant of time. Let t c be the delay of 
the ldf at to = o) c . The phase characteristic of the ldf is 

</>(w) = ~(t c + o3 c /a)o3 + 03 2 /(2a) + c, (28) 

where c is a constant, and a is defined in eq. (9). The ldf output is 

r ,,2 

y(t) = cos 



u c t + (t c + 03 c /a)(— co c ) + — - + c 

2a 



+ jfii™ 



03 c \ (C0 C + (i)/ 
(Oct + 03 jt — I T c H (tOc + 03 j) H h C 

a / 2a 



+ cos 



tO c < — W y Y — I T c H (We- 03 j) H 1- C 

a / 2a 



(29) 



In the above expression, various sidebands of the input am signal are 
delayed asymmetrically about the carrier frequency. This would, of 
course, distort the signal. After some tedious manipulations, the en- 
velope of y(t) is found to be 

Y(t) - | 1 + 2 mjcoso3j(t - Tc)cos(^- 



V mjCoso3j(t — Tc)sin( — ) 



2 \ -12-. 1/2 



(30) 



Comparing Y(t) and x(t), it is obvious that distortionless transmission 
results only if 03]/ (2a) «: 1. More importantly, the distortion in Y(t) is 
dependent on the input signal and, therefore, cannot be equalized 
easily. On the other hand, if synchronous demodulation is used in 
place of envelope detection, we obtain the first square bracket term in 
eq. (30), i.e., 



Y s (t) -1+2 77l,COSW, (t - Tc)C0s( -^- ) , 



03] 



(31) 



where the distortion shows up as cos[w y /(2a)] which is independent of 
the input signal, and equalization is thereby feasible. To do synchro- 
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nous demodulation, we need to know the instantaneous frequency of 
the carrier in y(t), the ldf output, which we have not addressed so far. 
Let us do so in the following. 
Stating eq. (16) again, 

y(t) = exp[y'27rO(0] A(t) exp 



T 2 



• exp[j2Tr(at - A/c - U + h)r\dr. (32) 

In the above we have, in part, an output chirp frequency of d®/dt = 
(U- at). The integral part, on the other hand, can be interpreted in 
two different ways: (i) It is a "quadratic" chirp transform of A{t). As 
such, little is known about this transform, (ii) It is the Fourier trans- 
form of A(t) exp[y'277-(/3 - a)T 2 /2], where A(t) can be viewed as an 
envelope modulation on a chirp signal with frequency (/? — ol)t. 
Equivalently, it is the convolution of the Fourier transform of A(t) 
and that of the chirp signal. The result may very well contain high 
frequency components depending on the magnitude of (/? — a) and the 
shape of A(t). In fact, simulation results show that the instantaneous 
frequency of y(t) can be quite different from (/ 4 - at). Therefore, 
frequency predictability is difficult in the general case, and synchro- 
nous detection as discussed above cannot be used easily. 

The following is a summary of the analytic results: 
(l) The mathematical expressions are complicated, and simulation 
is necessary to obtain numerical insights. 

(U) The compression operation can be justified by a quasi-station- 
ary model. Under this model and with envelope detection, multitone 
am leads to nonequalizable distortion. Furthermore, the quasi-station- 
ary model is valid only for small tcrs. 

(Hi) Synchronous detection is very difficult because of frequency 
unpredictability. 

III. SIMULATION 

This section describes simulation results. These numerical results 
show that the proposed technique indeed time compresses the input 
signal, but the output distortion is probably unacceptable for TV 
transmission using today's saw devices. 

3.1 Preliminary set-up 

A computer subroutine was written to simulate the time-compres- 
sion operation. It accepts an arbitrary input A (t) defined in the interval 
(0, Ti) and outputs y(t), y c (t), and y a (t), given by eqs. (15), (21), and 
(22), respectively. Note that y(t) is the rf output of the ldf, and y c {t) 
and y„(t) are the in-phase and quadrature components after synchro- 
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nous demodulation with <b(t) so that the ideal envelope detector 
output z(t) can be easily calculated from eq. (23). Also note that all 
three outputs have the same integral form, i.e., 

J = f(t)cos[2ir(at 2 + bt + c)]dt, (33) 

where we can assume t\ > 0, t 2 > t u and a > 0. In the computer 
program, we break the interval (0, T t ) into many small segments such 
that within each segment, a linear approximation of A{t) is valid. After 
doing so, the integration over each of these segments can be done in 
the form of eq. (33) with 

f(t)*mt + k. (34) 

The limits of integration t\ and t 2 in eq. (33) are, of course, the end 
points of the small time segment. Using the linear representation eq. 
(34) for f(t), I in eq. (33) can be integrated in closed form in terms of 
Fresnel integrals, and the result is 

/ - mlz + kh , (35) 

where 

7i = [ir/(2p)] 1/2 {[C(z 2 ) - C(zi)]cos B - [S(z 2 ) - S(zi)]sin B), (36) 

COS .o 

h = — — {sin xl - sin x 2 - q[v/(2p)] 1/2 [C(z 2 ) - C(2,)]} 
zp 

+ ^- {cos xl - cosx? + q[ir/(2p)] 1/2 [S(z 2 ) - S(2,)]}, (37) 



and 



p = 2ira; q = 2irb; r = 2ttc, (38) 



s = r -fe)' (39) 

Xi=>/pti+-2=, (40) 

2Vp 

x 2 = y/pt 2 + -^j=, (41) 

2Vp 

Xl 

«i— 7— . (42) 

vx/2 

x 2 
Z2 = -]=- (43) 
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The Fresnel integrals in eq. (36) are defined by 



'o 
and 



77 



C( Z )= I cos[^x 2 )dx (44) 



S(z)= | sm\^x 2 )dx. (45) 



There are also some simpler cases, e.g., a = or b = 0, which are not 
shown for brevity. Since time expansion can be achieved by reversing 
the delay slope of the ldf, the subroutine for simulating the time 
compression can also be used to simulate the time-expansion filter. 
Examples of both time compression and expansion will be shown later. 
In all subsequent simulations, the input is always 

Att |-/l. 0<<<T, (3l-Mf4 (46) 

A(l) [0, otherwise 

This rectangular pulse is, of course, not a representative of the video 
signal. However, except for edge "ringings," the system should com- 
press the pulse properly, and peak-to-peak (p-p) ripples at the center 
portion of the output pulse should give an indication of the magnitude 
of distortion involved. To examine the outputs y(t), y c (t), and y 8 (t), we 
have also developed a set of programs to estimate the instantaneous 
frequency at various time instants of y(t), as well as the p-p ripples of 
the output pulse. 

Because of the complexity of the equations and various computer 
routines, it is not easy to assure that there is no bug in the programs. 
However, we did make some runs for the special case a = /? (Fourier 
transform case) where results are known. The waveforms y c (t), y s (t), 
and z(t), and the output chirp frequency and offset frequency of the 
ldf given by eq. (31) were all verified carefully. Such a check assures 
the validity of the computer programs. 

3.2 Examples of time compression 

Four examples of time compression are discussed in this section. 
They all have tcr = 2 (i.e., 2:1 compression) and an input given by eq. 
(46) (see Fig. 5a). The chirp frequency range A/ c for the input x(t) to 
the ldf is (f 2 , /a), while the passband A/ of the ldf extends over 
( /i, f 4 ). The delays at U are denoted by t, (i = 1 to 4), respectively. The 
delay dispersion of the ldf is denned by At = n - t 4 . The descriptions 
for these examples are as follows: 

(i) The key parameters for the ldf are shown in Fig. 5b. In this 
example, we let f\ = U, fa = U and choose 
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Fig. 5 — Parameters for time-compression technique. Example 1 — (a) Input pulse (t ); 
(b) Delay of ldf; (c) Expected output (t). 



A/- c = A/"=60MHz. 



(47) 



From eq. (2), tcr = 2 with At = 32 /is. The time-bandwidth product 
(bt) of the ldf is defined by 



BT=(A/)(At). 



(48) 



Here, bt = 1920. Bandwidth product = 10,000 to 50,000 is considered 
as a practical range for saw filters. The "expected" output is illustrated 
in Figure 5c where the output compressed pulse appears between t = 
32 /is and 64 fis. The small ripples in the time intervals, to 32 /is and 
64 to 96 lis are to illustrate the nonzero output of the ldf in these 
regions. 

(ii) The ldf parameters are shown in Fig. 6a and the expected 
output in Fig. 6b. Here, A/" c = 60 MHz which is the same as in Example 
1 (Fig. 5), but At and A/" are both increased by a factor of 3. With 
A/= 180 MHz and At = 96 /is, bt = 17,280. 

(Hi) The parameters are shown in Fig. 7a, and the expected output 
is the same as that of Example 2 (Fig. 6) because the delay dispersion 
has remained the same. We increased A/ to 600 MHz yielding bt = 
57,600, which is probably a little beyond present-day state of the art. 
We used A/ c = 200 MHz. 

(iv) The parameters are shown in Figs. 6a and 6b. This is essentially 
the same as Example 3 in Fig. 7a, except A /"is increased to 1200 MHz, 
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Fig. 6 — Parameters for time-compression technique. Example 2 — (a) Delay of ldf; 
(b) Expected output (t). 

and bt = 115,200. This is probably not realizable in the near future. 
We used A/ c = 400 MHz. 

The above examples are arranged to illustrate the effect of increasing 
A/ (or the strengthening of the quasi-stationary model). 

With the input being a pulse 64-/xs long, the expected output for all 
examples is a compressed pulse, 32 /is in duration (tcr = 2). This is 
indeed so from the simulation results which show a compressed pulse 
approximately 32 ps long in the expected time slot. The output outside 
the compressed pulse duration is at least an order of magnitude lower. 
However, inside this compressed pulse, there are ripples created by 
the compression process, and these ripples are the resulting distortion 
that we are trying to estimate. The ripples are largest at the edges of 
the compressed pulse and smallest toward the center. Also, the ripples 
at the center are indications of the "best" performance of the time- 
compression filter. In our computer routines to estimate distortions, 
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we take some symmetric time interval, say ±AT d , about the center of 
the time-compressed pulse, and we record the maximum and minimum 
pulse magnitudes, denoted by 2 max and 2™, respectively. The p-p 
ripple distortion (over ±A7d) is defined by 



RD A 20 log 



- Zr 



0/2 



(49) 



\2max ■ ZuodI 

where rd is in dB. Both z max and Zmm are positive because A (t) is a 
positive pulse. We plot the p-p ripple distortion versus ATj in Fig. 8. 
The largest Aid is 16 /xs because that is the edge of the compressed 
pulse. Although it is difficult to translate the meaning of rd to a TV 
quality measure, it can be certain that a large rd (i.e., large ripple) 
would mean poor transmission quality. To make a TV system workable, 
an rd of less than —45 dB is probably necessary, although other 
applications may not require such a low distortion. 

Referring to Fig. 8, it is obvious that the distortion gets progressively 
smaller from examples 1 to 4. But the lowest distortions, despite the 
large bts involved, are still too excessive for high quality tv transmis- 
sion. Some more examples are provided in Appendix B for additional 
insight into the problem of ripple distortion. 

3.3 An example of time compression and expansion 

In this example, we use practical design parameters involving both 
time compression and time expansion (Fig. 9). Figure 9a shows a 64- 



-30 - 



40 




EXAMPLE 1 

-cr 



EXAMPLE 2 



2 4 6 8 10 12 14 16 18 

+ DURATION ABOUT CENTER OF OUTPUT PULSE IN MICROSECONDS 

Fig. 8 — Distortion results for time-compression examples. 
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Fig. 9— Parameters for (a) time-compression and (b) time-expansion examples. 

jus pulse as an input to the compressor, and the expected output. The 
actual ripple distortion (i.e., simulated result) over the 32-/xs com- 
pressed pulse is shown in Fig. 10a. To check the operation of the 
expander alone, we assume an input of a 32-jus pulse to the expander. 
The expected result is shown in Fig. 9b, and the simulated ripple 
distortion is plotted in Fig. 10b. It is clear from Figs. 10a and 10b that 
both the compressor and the expander lead to considerable distortion. 
To simulate a total system with compression and expansion, we 
used the distorted 32-jus pulse from the compressor as an input to the 
expander. The expanded output pulse is so distorted that it becomes 
meaningless to make a ripple distortion plot as before. Instead, we plot 
a small segment near the center of the expanded pulse in Fig. 10c to 
illustrate its excessive distortion. The performance of this total system 
is definitely not suitable for TV transmission. Some filterings were tried 
at the output of the expander, but little improvement was obtained. 

IV. CONCLUSION 

We have studied a time-compression technique based on a simple 
configuration of saw filters. The technique was derived heuristically 
and can be viewed as a quasi-stationary model for the chirp signals. 
Numerical results show that excessive distortion is created, and its 
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Fig. 10— Results for time-compression and time-expansion examples, (a) Ripple dis- 
"tion of time-compression example, (b) Ripple distortion of time-expansion example. 
Pulse ripple for time-compression and time-expansion examples. 

application to TV transmission is not suitable unless some kind of 
equalization is provided. One such equalization is the chirp transform 
processor 3 which involves considerably more complexity. Simpler 
equalizations may be possible but do not seem to be straightforward. 
As for other applications where the distortion requirement is less 
stringent, this approach may be feasible. 
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APPENDIX A 

Impulse Response of a Linear Dispersive Filter 

A linear dispersive filter (ldf) is a filter having a linear delay 
characteristic. Using the analytic-signal notation, its transfer function 
can be written as 

(Aco Aco 

i, -— So - at>£ - m 

0, elsewhere 
and 

A" Aco 

.D(co) = (t - scoo) + sco, — — < co - coo < — > ( 51 ) 

where G(co) and D(co) are the gain and group delay, respectively; to is 
the angular center frequency, and Aco is the bandwidth in radians; t is 
the delay at coo; and At is the delay dispersion over coo + Aco/2. The 
delay slope is defined by 

si^. (52) 

Aco 

Integrating D(co) with respect to co, we obtain the phase, 

s 9 Aco Aco 

4>(co) = -(to - scoo)co - - co , — — <co-coo<— , (od) 

where the integration constant has been dropped for convenience. 
Using exponential notation and neglecting all unimportant multiplying 
constants and constant phase shifts in subsequent discussions, the 
impulse response of the ldf is 



h(t) = Re I , exp 



Au 



Aw 



S 2 
(t — scoo)co + - CO 



exp(jut)du. (54) 



The above is simply an inverse Fourier transform of the transfer 
function. By a change of variable, it can be rewritten as 



h(t) = Re exp(y'coo£) . exp 



Au 



S 2 
-J | TOCO + - CO 



exp(jut)du. (55) 



There are two methods to solve for h(t) according to the above: (i) 
extend the limits of integration to ±oo and obtain a closed form solution 
easily; (ii) retain the finite integration limits and derive the result 
using Fresnel integrals. Both methods are shown briefly below. 
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Method I 



Changing the limits of integration to ±00 in eq. (55) and recognizing 
that the integration of the term exp(— j't co) leads to a delay of to in 
h(t), we obtain an integral in the form of 



-I 



exp ( —jkco 2 ) exp (j<i)t)du, 



(56) 



where k is a constant. Putting the integrand in Gaussian form by 
completing the square, we get 



J= \/-^exp 






Applying the above to eq. (55), we obtain the desired result 



or 



h{t) = Re exp(y<oo£)exp j 

(t - to) 2 



. (t - To) 5 



2s 



h(t) = cos 



coot + 



2s 



(57) 



(58) 



(59) 



However, because the original transfer function has a delay dispersion 
At defined about to, the valid range of t for eq. (59) is 



At At 

T0 __<,< T0 + _ 



(60) 



Method II 

Setting wi = wo — Aw/2 and Ti = to — At/2, eq. (55) is written as 



h(t) = Re exp(jcoit) I exp — j I tiw + - u 2 J 



exp(jcot)du). (61) 



Again recognizing that exp(-/Tia>) leads to a delay of n, we may 
substitute t = t — T\ and 

h(r) = Re expijuit) exp I -j - w 2 ) exp(jcot)du. (62) 

Completing the square in the integrand, we obtain 



&(t) = Re exp(y'wiT)exp [j — 



exp 



. s 

-7-1(0 

J 2 \ 2 



rfw. (63) 



The above can be integrated using Fresnel integrals, and the result is 
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h(r) = Reexp 1/ f «iT + g» J 

x{[C(y 2 ) - C(yi)] -j[S(y 2 ) - S(yi)]}, (64) 



where 

y 2 A 



-Awf 1-^-1. (M) 

7r V At. 



C(z) and S(z) are the Fresnel integrals defined in eqs. (44) and (45). 
Consider the bracketed term: 

[C(y 2 ) - C(yi)] -j[S(y 2 ) - S( yi )] A p(r)exp[je(r)] . (67) 

Neglecting small ripples, p(r) is constant over < t < At and vanishes 
outside this range. The phase term 0(t) is approximately constant over 
the same interval < t < At. Therefore, putting back the ti delay, we 
have 

'^- T1 ) + ( '- Tl)2 ' 



h(t) a Re exp \ j 



2s 



(68) 



APPENDIX B 

Additional Examples On Ripple Distortion 

Additional examples are provided here for more insight into the 
phenomenon of ripple distortion. The first case of interest is that of 
band-limiting effect on the input pulse. We use a 32-/XS input pulse and 
low-pass filter it with a raised-cosine characteristic, where the gain is 
unity from zero to 8 MHz, 0.5 at 9 MHz, zero at 10 MHz, and the delay 
is constant over the passband. The output is, of course, a 32-/xs pulse 
with ripples. The magnitude of these ripples are plotted in Fig. 11 in 
the manner similar to Figure 8. 

The second case of interest is that of ripple distortion caused by the 
linear delay slope of the ldf. Three specific examples are provided to 
illustrate this effect: 
Case A: The input to the time-compression filter is a 32-jus pulse. It 

is modulated by a CW frequency of fo = 1600 MHz. The ldf 

has a delay characteristic as shown in Fig. 1, where fo = 1600 

MHz and A/= 2400 MHz. 
Case B: The conditions are identical to those of Case A, except fo = 

1000 MHz and Af= 1200 MHz. 
Case C: The conditions are identical to those of Case A, except f = 

700 MHz and A/ = 600 MHz. 
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Fig. 11 — Distortions because of bandlimiting and delay slope. 

All three examples above illustrate the case of no time compression, 
but simply a constant delay through the ldf as viewed by the quasi- 
stationary model. Note that the magnitude of the delay slope increases 
by a factor of two from each example to the next and, hence, an 
increase in ripple distortion as seen from the results plotted in Fig. 11. 
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